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ABSTRACT 

We explore how the co-evolution of massive black holes (MBHs) and galaxies is affected by environ- 
mental effects, addressing in particular MBHs hosted in the central galaxies of clusters (we will refer to 
these galaxies in general as 'CGs'). Recently the sample of MBHs in CGs with dynamically measured 
masses has increased, and it has been suggested that these MBH masses (Mbh) deviate from the 
expected correlations with velocity dispersion (a) and mass of the bulge (Mb u i ge ) of the host galaxy: 
MBHs in CGs appear to be 'over-massive'. This discrepancy is more pronounced when considering the 
Mbh — cr relation than the Mbh — Mb u i ge one. We show that this behavior stems from a combination 
of two natural factors, (i) that CGs experience more mergers involving spheroidal galaxies and their 
MBHs, and (ii) that such mergers are preferentially gas-poor. We use a combination of analytical and 
semi-analytical models to investigate the MBH-galaxy co-evolution in different environments and find 
that the combination of these two factors explains the trends observed in current data-sets. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: formation 
— black hole physics 



1. INTRODUCTION 

The discovery of correlations between MBHs and their 
hosts flMagorrian et all 119981 iFerrarese fc Merritt] 12001 
IGebhardt eHdlE pOO) has been taken as one of the main 
elements in support of a co-evolution between MBHs and 
galaxies, that in turn brought to the suggestion that en- 
ergy input from an accreting MBHs, as quasar or Ac- 
tive Galactic Nucleus (AGN) may regulate star forma- 
tion in the host, or create a symbiosis between MBH 
and stellar growth. Indeed, theoretically there is rea- 
son to expect that energy or momentum driven out- 
flows establish correlati ons that scale as Mbh oc cr° an d 
M BH oc ct 4 respectively (|Silk fc Reed[T^lFabian|[l999l ). 
The correlations between MBHs and g alaxies, derived 
on a sample of about 50-70 MBHs (e.g., Gulteki n et al.1 
[20091: iGraham et all [201lt IMcConnell fe Mai 12012ft "are 
also extrapolated to the whole population, assuming that 
each galaxy hosts a MBH consistent with the correlation 
and its scatter to derive global pr operties, such as th e 
mass density in MBHs today (e.g.. iShankar et aLH "2004). 
These correlations, or deviations from such correlations, 
also appear to hold the key to understandin g the for- 
matio n and evolution of the MBH population (|Volonterl 
2012, and references therein). 

Recent measurements of the masses of MBHs in CGs 
(we include in this definition central dominant galaxies 
as well as brightest cluster galaxies) seem to suggest that 
these MBHs arc 'over-massive' compared to expectations 
from of the Mbh — cr relation, but they appear more con- 
sistent , jdbejtwithji^^ 



lation (jMcConnell et a ll 12011 . 
(|20T1 and IGraham fc Scott 



IMcConnell fc Mai 
fit for the Mbh — cr 



correlation find that the relation becomes much steeper, 
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e.g, from M BH oc cr 4 2 ([Giiltekin et al.ll2009D to M BH oc 
cr 5 - 6 . The same does not happen for the Mbh — Mb u i ge re- 
lation, that appears to be consistent with being linear in 
both cases (although including MBHS in CGs increases 
the normalization). There is also tentative evidence for 
over-massive MBHs in CGs from comparisons with the 
expectati ons from the Fundamental Pla ne of black hole 
activity (Hlavacck- Larrondo et all I2012T ) . Feedback af- 
fectin g galaxies in different w ays depending on potential 
well (Boo th fc Schavdl2011ft an d environment, via gas- 
rich mergers and cooling flows (Zubovas & King 120121) 
has been proposed to be partially responsible for the 
properties of these MBHs. 

We investigate here two guiding ideas related to the 
influence of mergers between spheroidal, gas-poor galax- 
ies ('dry' mergers) and of MBH-MBH mergers to de- 
termine the astrophysical drivers of the possibly differ- 
ent relationship between MBHs and their hosts in the 
case of CGs. It is well established from the theoretical 
point of view that (parabolic) dry mergers consistently 
grow a galaxy's mass, luminosity and radius more than 
a gal a xy's velocity dispe r sion (e.g.. ICiotti fc van Albada 



2001 



2009; 



Ciotti et all 120071 iNaab et all 120091: iNipoti et al 



Shankar et al.l 120111: IQser et all 120121: IHilz e~ 



20121 ) . especially w hen a galaxy is dominant over th e 



general population ([Ciotti et al.ll2007l : |Naab_eLaLj[2009) . 
Qualitatively this implies that eventually the largest 
galaxies deviate from the Mbh — cr relation. Therefore, 
if CGs grew predom inantl y through dry mergers (e.g . , 
lOstriker fc Tremainel Il975l lHausman fc Ostrikerl [19781) . 
and their MBHs grew mostly through MB H-MBH merg- 
ers ([Malbon et al.l 120071 : lYoo et al.l 120071 ) in such dry 
mergers then the expectation is th at they will be outliers 
in th e Mbh — cr re l ation ( see also iBovlan-Kolchin et al.l 
[2001 lLauer et al.l [20071 iZhang et al.l I2012D . Sup- 
port for the influence of dry mergers in shaping the 
structural properties of CGs comes from their steeper 
radius-luminos ity relation and flatter Faber- Jackson 
relation (e.g., lLauer et afl 120071 iBernardi et al.l 120071 : 
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iDesroches et~aTl [20071 ) . Additionally, if MBHs in CGs 
experience a significant mass increase because of MBH- 
MBH mergers with respect to MBHs in a field environ- 
ment, the extra boost in mass may also contribute to 
explain the properties of MBHs in CGs, if gas accretion 
is responsible for establishing the relationships in the first 
place. 

We present here models of MBH growth in galaxies 
that take into account both the cosmic environment, i.e., 
the frequency and properties of galaxy and MBH merg- 
ers, as well as the evolution of the virial properties of 
galaxies during galaxy mergers, with particular attention 
to mergers between gas-poor spheroids. 

2. MERGERS AND BLACK HOLE GROWTH 

We us e the latest incarnat ion of our semi-analytical 
models ([Volonteri et al.|[2012f) to provide a simple esti- 
mate of the influence of dry mergers on the build-up of 
MBHs and their hosts in different environments. We 
summarize here the main assumptions. We keep our 
model for galaxy morphology and gas content as sim- 
ple as possible. Morphology is related to the merger 
history, using a three-parameter model, where spheroid 
(we equivalently refer to spheroids and 'gas-poor' galax- 
ies in the following) formation depends on both halo 
mass ratio and the absolute halo mass, and a spheroid 
can re-acquire a disc th rough cold fl ows a nd mergers 
with gas-rich galaxies. iKoda et "all (|2009| ) show that 
the fraction of disc- vs spheroid-dominated galaxies is 
well explained if the only merger events that lead to 
spheroid formation have mass ratio >0.3 and virial ve- 
locity > SSk ms" 1 ; also, the merger tim escale (calcu- 
lated following iBovlan-Kolchin et al.l 120081 ) must be less 
than the time between when the merger starts and today, 
z = 0. We assume that spheroids form after a merger 
that meets these requirement. We additionally allow gas 
to re-condense ('gas-rich' galaxy) after 5 Gyrs in galaxies 
with virial velocity < 300 km s^ 1 where no major merg- 
ers occurred to include the effect of gas accretion and 
cold flows. 

We wish to keep our models as simple as possible, 
while making sure that the properties of the MBHs we 
study are correctly determined through the cosmic evo- 
lution of their hosts. We do not explicitly model the 
evolution of the baryonic component of the host galax- 
ies throug h cooling, star format i on and various feed- 
backs (see iFanidakis et all 120 lit iFontanot et al.l 120111 : 
iHirschmann et al.ll2012l and references therein for mod- 
els that treat in detail semi-analytically the baryonic 
component of galaxies and its link to MBH evolution). In 
our models we use only one parameter to link the host 
halo to the central MBH, and it is the halo's velocity 
dispersion. We link the velocity dispersion of the halo to 
the asymptotic virial velocity (i> V ir) assuming a spherical, 
isothermal halo, so that cr vir = u vir /\/2. We calculate the 
circular velocity from the mass of the host halo and its 
redshift. 

At high redshift we seed dark matter halos with MBHs 
created by gas collapse. Specifically, we adopt here 
the formation model detailed in Natarajan & Volonteri 
(2011) based on Toomre instabilities (Lodato & Natara- 
jan 2006). 

Our model includes MBH-MBH mergers, merger- 



Table 1 



M halo 


Dry mergers 


Dry mergers 


A^BH,o/-^acc 


(M e ) 


with MBHs 


(total) 




10 15 


4.2 (cr n =2.0) 


4.6 ((Tn=2.1) 


10.9 ± 9.7 


10 14 


1.7 (cr„=1.3) 


1.8 K=1.4) 


2.4 ± 1.0 


4 x 10 13 


0.3 (cr„=0.5) 


0.4 (<r n =0.6) 


1.9 ± 0.7 


2 x 10 13 


0.2 (<r„=0.4) 


0.3 (cr„=0.5) 


2.0 ± 0.7 


10 13 


0.1 (cr„=0.2) 


0.1 (<r n =0.2) 


1.8 ± 1.8 



Note. — Number of dry mergers (including their 1 — a variance) 
experienced by the central galaxy, regardless of its morpholgy, of 
a halo that by z = has the mass listed in Column 1, and growth 
channels for its MBH. 



driven gas accretion, stochastic fueling of MBHs through 
molecular cloud capture, and a basic implementation 
of accretion o f recyc led gas. Details are given in 
iVolonteri et al.l (|2012l) . Despite its simplicity, our ap- 
proach produces a population of MBHs and AGN con- 
sistent with the observed one in terms of, e.g., luminosity 
function of AGN, relationship between MBHs and their 
hosts, high-redshift quasars. 

In Table 1 we list the number of dry mergers experi- 
enced over its lifetime by the central galaxy of a halo, as 
a function of the halo mass at z = (the average is per- 
formed over 20 to 35 different halos for each halo mass). 
We can consider the 10 15 M halo as a cluster-sized halo, 
the 10 14 M halo as a group-sized halo, and the smaller 
halos as galaxy-sized ones. The dry-ness of a galaxy dur- 
ing a merger is determined in a very simple way, i.e. we 
require that both merging galaxies are gas-poor. 

As expected, the number of dry mergers decreases as 
the halo mass decreases, therefore their influence will be 
strongest for central galaxies in clusters, and possibly in 
groups. In Fig. [T] we focus on the MBH growth in six 
different realizations of 10 15 M halos. For instance, the 
growth of the two MBHs in the top panels is strikingly 
different. In one case it is dominated by MBH-MBH 
mergers at z < 2, while the second MBH does not expe- 
rience any important MBH-MBH merger at late cosmic 
time. If we define the relative growth through mergers as 
the difference between the MBH mass at z = 0, Mbh,o, 
and the cumulative accreted mass, M acc , we can esti- 
mate the importance of MBH-MBH mergers in a MBH's 
history. In Table 1 we report this information in the 
fourth column. The definition we adopt cannot account 
for accretion on the MBHs that merge with the one 
in the central galaxy. Therefore, the fact that we find 
A^BH.o/Afacc > 1 does not mean that MBHs grow pre- 
dominantly through mergers, and in fact if we removed 
accretion from our models the final MBH mass would be 
about two orders of magnitude smaller. Overall, there- 
fore, MBHs grow through accretion o f gas, in a way 
consi stent with independent estimates (|Yu fc Tremaind 
[2001 . 

To test the robustness of our results we changed the 
MBH feeding scheme, by assuming that all MBH accre- 
tion activity is driven by galaxy mergers, that the accre- 
tion rate is fixed to 30% of the Eddington rate and that 
accretion stops once the MBHs have reached the value of 

Mbh = 10 8 (ervir/aOOkms- 1 ) 4 M . We find that while 
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Figure 1. Examples of MBH growth in six randomly chosen cen- 
tral galaxies of 10 f5 M halos. The red solid curve shows the MBH 
mass at a given time, the dark green dashed curve shows the cumu- 
lative mass gained in accretion events, blue asterisks mark MBH- 
MBH mergers with mass ratio > 1 : 10. As an example, the growth 
of the MBH in the top-left panel is dominated by MBH-MBH merg- 
ers, while that of the MBH in the top-right panel is dominated by 
accretion through various channels (merger driven in gas-rich merg- 
ers, through accretion of gas clouds, and through recycled gas). 



small quantitative differences exist, qualitatively the re- 
sults of our investigation are unchanged. We conclude 
that, if the Mbh — & relation is established because of 
accre tion-driven feedback (e.g-. lSilk fc Re es 1998; Fabian 
1999), MBHs in CGs may not necessarily obey the same 
the Mbh ~ G relation, as their masses may grow sub- 
stantially through mergers after accretion dwindles. 

3. MERGERS AND GALAXY STRUCTURE 

From the evolutionary histories of our models we ex- 
tract the series of mergers that the central galaxy of 
the main halo experiences from z = 1. We record the 
MBH mass in the merging galaxies and the dark mat- 
ter halo masses, M^.i and (here onwards the sub- 
scripts refer to each one of the two galaxies, with the 
convention that subscript 1 labels the central galaxy of 
the main halo, which is not necessarily the most mas- 
sive galaxy of the merging pair). To each halo at the 
starting redshift we assign a stellar mass, M*, using 
the fits by Behroozi et al. (2010, see also Nipoti ct 
al. 2012), a projected velocity dispersion velocity dis- 
persion, cr, consistent with the Faber- Jackson relation, 
and an effectiv e radius determined throu gh the Funda- 
mental Plane (jBinnev fc Tremaind [20081 pages 23-24), 
including scatter in all relations. We included a redshift 
dependence of the Faber- Jackson an d Kormendy rela - 
tions, using the scalings suggested bv lOser et al.l (|2012t ). 



namely that a oc (1 + z)° and R c oc (1 + z) (see 
also Nipoti et al. 2012). We then calculate the M* and a 
resulting from the merger following Ciotti et al. (2007) 
as follows. We account for weak homology by relating cr 



to the virial velocity dispersion, 

24.31 + 1.91n + n 2 



a v 44.23 + 0.025n + 0.99n 2 ' 
r v ^ 250.26 + 7.15n 
FL 



77.73 + n 2 



(2) 



adequate for Sersic index n, 2<n<12. We assume 
that the Sersic index of the resulting galaxy is n = 
1 + max(ni,ii2), where n\ and n% are the Sersic in- 
dices of the progenitors (each galaxy starts with n deter- 
mined from the virial coefficient, K v (n) = GM*j 'R c a 2 — 
(r v /R c ) x (a v /a) 2 ). 

In case of mergers between gas- rich galaxies or between 
a gas-poor and a gas-rich galaxy, we assume that each 
gas-rich galaxy has a gas mass M g = aM*, with a = 4 
and that a fraction rj = 0.05 (My/M/ji) of the gas is 
converted into stars: 



M* = M*i + A/*2 + Tf(M gl + M g2 ). 



(3) 



We find the velocity dispersion of the newly formed 
galaxjQ as: 



a:. 



Mtotl „ 2 , M tot2 , 

-Aoo. 



Aft, 



(4) 



where M totl = Af*i + M gl , M tot2 = Af* 2 + M g2 , M tot 
Aftoti + Af to t2, and 



= 1 



1 + ai 



(5) 



and a similar expression holds for A 2 . From the second 
merger of the sequence onwards, the central galaxy re- 
tains the properties derived in the previous step, while we 
assign M*2 and cr„2 to the merging galaxy as described 
above. Throughout our experiment, the MBH mass is de- 
termined from the semi-analytical model, that includes 
both accretion and MBH-MBH mergers. We note that 
not all galaxies host MBHs, i.e., the merging galaxy may 
or may not contribute to the MBH growth. 

We present the results of this experiment in Fig. [5J 
In the top two panels we focus on MBHs hosted in the 
central galaxy of 1O 15 M halos at z — (akin to CGs), 
while in the bottom panel the focus is the central galaxy 
of a galaxy-sized halo (1O 13 M at z = 0). We assessed 
that the results at z = are qualitatively unchanged if 
we ignore the redshift evolution in the Faber- Jackson and 
Kormendy relations, with MBHs and galaxies occupying 
the same region in the Mbh — & relation. CGs are char- 
acterized by MBHs that consistently deviate from the 
expected correlations, being over-massive at fixed galaxy 
properties, occupying the same range as observations. 
MBHs in elliptical galaxies that are not CGs, where the 
influence of MBH-MBH mergers and dry mergers is much 
milder tend instead to sit closer to the global MBH-host 
correlations. This result can be interpreted as overall 
steeper and higher normalized relations with respect to 
previous estimates (compare solid and dashed black lines 
in Fig. 2). At z < 1 most CGs are the dominant galaxies, 

4 The presence of dark matter can be included through an extra 
Cdm parameter in our scheme (see Ciotti et al. 2007), assuming 
that stars and dark matter are similarly distributed. Inclusion of 
the dark matter halo represents a small correction as long as we 
are interested in the region within the effective radius. 
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Figure 2. Examples of Mbh — & and Mbh — A^bulge in 20 ran- 
domly chosen central galaxies of 10 15 M Q (top) and 10 13 Mq (bot- 
tom ) halos. The gray erro rbars are observed MBHs and galax- 
ies (McConnell & Ma 2012, in the top pan el we show CGs only) , 
the solid blue lines a r e the fits derived by Giiltckin ct al.. (120091) 
and IMa rconi fc Huntl 12003 ). the dashed lines are the fits from 
IMcConnell fc Mai l|2012l ). The red dots show the location of model 
MBHs at z = 0. The dark red curves show two examples from the 
evolutionary histories: horizontal rightward swings in the Mgu—a 
panels occur when the galaxy merges with a more massive galaxy, 
vertical (leftward) movement is characteristic of dry mergers with 
similar (smaller) galaxies. 



and they merge with smaller galaxies that consistently 
have o V 2 < cr v i and therefore their velocity dispersion 
cannot increase (see Eq. 2|) . On the other hand a non-CG 
galaxy has a higher chance of merging with a galaxy with 
a V 2 > (Jvi, with the merger remnant having a v > <j v \ (see 
the dark red tracks in Fig. 2 for an ex ample). 

We note (see also lNipoti et al.ll2012[ ) that this scheme 
tends to overproduce stellar masses by z = 0. In fact the 
M star — Mh relationship peaks at = 10 12 M®, there- 
fore by merging galaxies close to the peak the remnant 
galaxy ends up having an increased M sta r — Mh- While 
our merger sequence includes both galaxies below, at and 
above the peak, we find that we consistently overpredict 
stellar masses at z — with respect to the scaling we 
would obtain directly from the M sta r — M^ relationship, 
and that in general bulges appear too massive (this is evi- 
dent in the right panels of Fig. [2]). This problem would be 
alleviated if we i ncluded corrections for mass lost in the 
merging process (jNipoti et al.l [20 03). We also consider 
parabolic mergers only, where energy is perfectly con- 
served. However, not all mergers involving a CG are nec- 
essarily parabolic. On the one hand, since the galaxies 
move in the potential of the cl uster, hyperbolic me rgers 
may occur. On the other hand, iNipoti et al.l (|2003l ) note 
the effect of dynamical friction, that braking the galaxy's 
orbit may induce elliptical merging. Hyperbolic and el- 
liptical mergers have competitive effects on the evolution 
of ct. Elliptical mergers with negative orbital energy in- 
crease the final ct, and viceversa. This effect, as well as 



the suggestion by iHilz et ail ([20121 ) that energy transfer 
from bulge to halo grows the velocity dispersion, would 
improve the match of our exercise with observations. 

4. CONCLUSIONS 

In this Letter we have highlighted the effects that merg- 
ers have on the MBH population in CGs. Two main 
factors contribute to their evolution. Firstly, CGs expe- 
rience many more dry mergers with spheroids. Parabolic 
dry mergers grow a galaxy's mass, luminosit y and radius 
more than a galaxy velocity dispersion (e.g.. Ciotti et alJ 
2007 [ INipoti et al.ll2009t IShankar et aDboili: IOser"et~aTI 
2012t Iffilz et al.ll2012D . If in a given merger M BH and 



A^buigc increase relatively more than ct (see the discus- 
sion in Ciotti et al. 2007), a sequence of such mergers 
will lead to more massive MBHs at fixed ct. In Fig. 2 
(left panels) this corresponds to moving upward faster 
than rightward. 

Secondly, the sheer number and mass contribution 
of MBH-MBH mergers occurring in CGs galaxies is 
much higher (see Table 1). If we assume that correla- 
tions between MBHs a nd hosts are estab l ished because 
of AGN feedback (e g ISilk fc Reesl [T99l iFabianl fiMk 
IDi Matteo eTall200a Hopkins et aLlhooa ). then if MBH 
mergers contribute to the MBH growth after the bulk 
of quasar/AGN activity has ceased (Fig. 1), the MBH 
mass increase brought by these mergers will push the 
MBH upwards and out from the correlation established 
through feedback. One important caveat, however, is 
whether MBH binaries can merge efficiently in gas-poor 
environme nts, because of the so-c alled 'final parsec prob- 
lem' (e.g.. iBegelman et al.lll986T ), although various ef- 
fects, such as triaxiality and rotation, as well as the pres- 
ence of mass ive perturbers. may i ncrease the orbital de- 
cay rate, see iColpi fe Dottil ([20111 ) for a recent review. 

Our models are at variance with other models that 
study the impact of M BH mergers on the est ablish- 
ment of correlations (e.g., Jahnke & Maccio 2011) as we 
do not assume that MBHs populate all galaxies. In 
fact, the presence or absence of a central MBH leads 
to different evolutionary paths in the Mbh ~ c and 
Mbh — -^buige relations. We can consider two extreme 
cases. Let us assume that all galaxies host a MBH. Then 
at each merger M* and Mbh increase as the sum of 
those in t he two galaxies, barring for the effect s of stellar 
escapers ([Nipoti et all [20031 ; IHilz et al.l [Ml and non- 
linea r addition of MB H masses (|Ciotti fc van Albadal 
120011: ICiotti et al.l 120071 ). However, Eq. 3 shows that ct 
would stay the same or slightly decrease. This corre- 
sponds to a vertical upward movement in the Mbh — & 
plot, eventually leading to MBHs that are over-massive 
for their ct, but are not outliers in the Mbh — -^buige cor- 
relation. The other extreme case assumes that only the 
main galaxy hosts a MBH. Then at each merger M* in- 
creases, while Mbh and ct do not. Eventually the MBH 
in the galaxy that results from the merger sequence will 
be under-massive for its bulge mass, but it will not be 
an outlier in the Mbh — & relation. 

Broadly speaking, the models of MBH evolution that 
we adopt for this Letter (Volonteri et al. 2012) predict 
that the most massive MBHs, except for those hosted 
in CG s, are those that are best c orrelated with their 
hosts ([Volonteri fc N ataraian| [2009|) if their build-up is 
driven by a combination of accretion and mergers that 
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includes both gas-rich and gas-poor galaxies, and dis- 
persion should increa se at low MBH/galaxy masses (see 
Fig. 2, top panel in iVolonteri et al.l [20121) . where the 
MBH mass, even at z = tr aces the properties of 
the M BH formation mechanism (Ivan Wasse nhove et all 
[201(1 . It may well be that if different processes shape 
the MBH mass at different galaxy masses (MBH forma- 
tion at the lowest masses, AGN feedback at intermediate 
masses, MBH and dry mergers at the highest masses), 
there is not a unique link that straddles throughout the 
whole range. 
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